Tkachenko Oscillations and the Compressibility of a Rotating Bose Gas 
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The elastic oscillations of the vortex lattice of a cold Bose gas (Tkachenko modes) are shown to play 
a crucial role in the saturation of the compressibility sum rule, as a consequence of the hybridization 
with the longitudinal degrees of freedom. The presence of the vortex lattice is responsible for a q 2 
behavior of the static structure factor at small wavevectors q, which implies the absence of long 
range order in 2D configurations at zero temperature. Sum rules are used to calculate the Tkachenko 
frequency in the presence of harmonic trapping. Results are derived in the Thomas-Fermi regime 
and compared with experiments as well as with previous theoretical estimates. 
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In 1966 Tkachenko developed the theory of the elas- 
tic oscillations of a vortical lattice in incompressible su- 
perfluids, predicting the dispersion law lot = y/Ml/imq 
where m is the mass of the fluid particles, q is the 
wavevector of the wave and Q is the angular velocity of 
the fluid, related to the number n v of vortices per unit 
surface by the relation n v = 2£lm/h. The frequency lot 
is calculated in the reference frame rotating with angu- 
lar velocity Q. The Tkachenko modes are peculiar of 
superfluids, where the formation of singly quantized vor- 
tices in the rotating liquid gives rise to regular crystalline 
structures (Abrikosov lattice) , stationary in the rotating 
frame. The theory of the Tkachenko modes was later de- 
veloped by Baym Q , who generalized the hydrodynamic 
theory of superfluids taking into account the elasticity of 
the vortex lattice. The effects of the compressibility on 
the Tkachenko modes were discussed in details by Sonin 
The compressibility changes dramatically the dis- 
persion law at small wavevectors q. In particular, when 
q <C fi/c, where c is the sound velocity, the dispersion 
relation is no longer linear, but becomes quadratic in q 
as a result of the hybridization with the sound waves. 

The availability of vortex lattices in rapidly rotating 
Bose-Einstein condensates 0, has stimulated new the- 
oretical studies 0,0 of the Tkachenko oscillations in har- 
monically trapped atomic gases. In particular in Ref. 6] 
the discretized Tkachenko frequencies have been calcu- 
lated taking into account the finite and inhomogeneous 
nature of the system. In Ref. Q the corresponding values 
of q have been employed within an improved dispersion 
law which includes the effects of compressibility. Fully 
numerical calculations of the Tkachencko waves in har- 
monically trapped gases have been also recently carried 
out using Gross-Pitaevskii theory ||. The first experi- 
mental observation of the Tkachenko modes in harmon- 
ically trapped gases has been recently reported by the 
group of JILA 9] . These modes represent the low energy 
counterpart of the hydrodynamic modes already investi- 
gated experimentally in the presence of the vortex lattice 
|3, hi good agreement with theory [r(ij . 



The main purpose of the present work is to show that, 
despite their elastic nature, the Tkachenko oscillations 
can be naturally excited using density perturbations, 
thereby opening new perspectives of experimental inves- 
tigation. In particular we will show that, in the homoge- 
neous case, these modes exhaust the compressibility sum 
rule in the limit of long wavelengths. In harmonically 
trapped configurations the sum rule approach will be 
used to calculate the frequency of the lowest azimuthally 
symmetric Tkachenko mode and explicit results will be 
derived in the Thomas-Fermi regime. 

Let us start our discussion by evaluating the density 
response function of a uniform gas containing a vortex 
lattice rotating in the x-y plane with angular velocity 
f2 = Q e z at zero temperature, where e z is the unit vector 
along the z-direction. In a compressible fluid the condi- 
tion of uniformity can be fulfilled by adding the harmonic 
potential term mfl 2 (x 2 + y 2 )/2, which compensates the 
centrifugal effect produced by the rotation. This corre- 
sponds, in the rotating frame, to using the Hamiltonian 
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where N is the number of atoms, p is the momentum 
and g is the coupling constant of the 2-body interaction. 
The density response is easily evaluated using the coupled 
hydrodynamic-elastic formalism 0, 01 characterized, in 
addition to the hydrodynamic energy functional 
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where v is the velocity field in the laboratory frame and 
n is the density, by the elastic term 
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sensitive to the deformation of the lattice through the 
vortex displacement field e and characterized by the 
clastic parameters C\ and C2. At equilibrium one has 
v = f2 A r and e = 0. The dynamic response function 
takes the form 
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with r\ — > + , where oj+, u>- are, respectively, the upper 
and the lower branches of the energy spectrum. Expres- 
sion Q holds in the macroscopic regime q <C (H/mfl)^ 1 ^ 2 
corresponding to wavelenghts larger than the average dis- 
tance between vortices. The general expression for lu± as 
a function of C\ and Ci has been derived in Ref. Q. 
Here we report the results in the Thomas-Fermi regime 
Ml <C mc 2 , corresponding to the condition that the size 
h/mc of the vortex cores is much smaller than y/ h/rMl. 
In the Thomas- Fermi regime one has C2 = — C\ = nMl/8 
7J. In this case the upper branch follows the disper- 
sion law u} 2 ^ = 4S1 2 + c 2 q 2 and exhibits a gap at q = 0. 
Conversely the low frequency branch, hereafter called 
Tkachenko branch (w_ = lot), obeys the gapless law 
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For large q Eq. (J5J reproduces the original Tkachenko law 
y/ Ml /4m (7, while for small g it exhibits the quadratic be- 
havior lot = \J h / IQmOu cq 2 . The transition between the 
q 2 and q dependence takes place at values q ~ O/c which, 
in trapped condensates, can be significantly larger than 
the inverse of the radial size of the system. This suggests 
that the effects of compressibility, characterizing the q 2 
dependence, play a crucial role in the Tkachenko modes 
of a trapped gas, as explicitly pointed out in Ref. 0. 

Starting from the density response function J3J one 
can easily calculate the energy weighted and the in- 
verse energy weighted sum rules relative to the density 
operator p q — J2k=i e ~ lqXk ■ These are given, respec- 



tively, by mi(pg) = J2 n \( n \Pq\ty\ 2E nO and m_i (/£>,) 



J2 n \( n \Pq\^)\ 2 ^n0y where J2 n is the sum over all the ex- 
cited states and E n0 = E n — E is the difference between 
the eigenenergies of the excited state \n) and of the ini- 
tial configuration |0) containing the vortex lattice. The 
relation between these sum rules and the asymptotic be- 
havior of the density response function is given by [l"H 
x(q,u)uj^oo = -2mi(p q ) I \Mj) 2 and x(g,0) = 2rn-i(p q ). 
In the first case one recovers the model independent /- 
sum rule m\(p q ) — N(hq) 2 /2m which, at small q, is 
exhausted by the high energy branch lo + . By taking 
the lu — > limit of x(?>w) one instead finds the result 
m-i(pq) = N/2mc 2 also known as the compressibility 
sum rule. The Tkachenko branch plays a crucial role in 
satisfying the latter sum rule. In fact, because of the gap, 
the high energy branch lo + contributes to m-i(p q ) only 
through terms of order q 2 . 



It is also worth noticing that the static structure factor 
S(q) = N^ 1 J2 n \ { n \Pq\®}\ 2 is deeply affected by the rota- 
tion of the gas and behaves like q 2 , differently from what 
happens in non rotating interacting fluids where it is lin- 
ear in q \l2\ . By using the zero temperature relationship 
NS(q) — (ft/ 71 ") Jo°° I m x( w ) dw one finds the result 
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when q ~ * 0, where the second term in the parenthesis 
is the contribution of the Tkachenko branch 13] . The 
corresponding suppression of the density fluctuations re- 
sults in a dramatic enhancement of the fluctuations of the 
phase of the order parameter, which destroy long range 
order. This can be easily seen using the uncertainty prin- 
ciple inequality 2S(q)(2n q + 1) > no [l3 |. where n q is the 
particle occupation number and no is the Bose-Einstein 
condensate fraction. Since S(q) — > q 2 one finds that n q 
diverges at least like 1/g 2 at low q, thereby ruling out 
Bose-Einstein condensation in 2D even at zero tempera- 
ture [iH . 

In the following we will use the sum rule technique to 
evaluate the Tkachenko frequency through the ratio 
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between the energy weighted and inverse energy weighted 
sum rules relative to the excitation operator F. Our final 
goal is to derive explicit results in the presence of har- 
monic trapping. Eq. (JJJ provides a rigorous upper bound 
to the frequency of the lowest energy mode excited by F. 
The proper choice of the operator is a crucial step in the 
calculation. For example, in the uniform case it would 
not be appropriate to use the density operator p q for F 
since, as already pointed out, the /-sum rule mi{p q ) is 
exhausted by the upper branch at small q and the ratio 
would not coincide with the Tkachenko frequency. 

The general strategy for the identification of the excita- 
tion operator is suggested by the fact that the Tkachenko 
modes have zero energy in the hydrodynamic approx- 
imation ^(|, where elasticity effects are ignored. So 
one should look for excitation operators F whose en- 
ergy weighted sum rule mi(F) = ([F* , [H, F]})/2 vanishes 
when evaluated in the hydrodynamic approximation. 

In a uniform system this condition is satisfied by the 
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mClx)]k- In fact the corresponding double commutator 
takes the form [F*,[H,F]] - (ft 2 g 4 /4m 3 2 ) (p„ - 

infix) 2 , and its expectation value identically vanishes if 
one uses the hydrodynamic prescription p = mv[r) with 
the equilibrium condition v = f2 A r. The elastic contri- 
bution to mi (F) can be conveniently calculated applying 
to the equilibrium configuration the unitary transforma- 
tion U = e l9S , where S = (F + F^)/2 and 9 is a small pa- 
rameter. Due to the presence of the non local transverse 
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term e~ tqx py in F, the transformation U gives rise to 
the vortex displacement e(9) — — 6(Hq/2mSl) s'm(qx)e y , 
where e y is the unit vector along the y-direction. By 
calculating the corresponding clastic energy change @ 
and considering the expansion E{&) = (U^^-HU) ~ 
E Q + 9 2 ([S, [H, S]])/2 for the total energy of the system, 
one derives the result m\{F) = N(hfl/4:m)(h 2 q 4 /8f2 2 m). 

To obtain m-i(F) we calculate the static response 
x(F) to the perturbation — XF + h.c. using the hydro- 
dynamic energy functional and the relationship x(F) = 
2m_i(-F l ). The elastic term provides higher order cor- 
rections. For small values of q the leading contribution 
arises from the density component p q in the Tkachenko 
operator and is fixed by the compressibility of the gas. 
For larger values of q, or for incompressible fluids, the 
static response is instead determined by the transverse 
current term e~ lqx (p y — mClx). In general one finds the 
result to-i(F) = (N /%VL 2 mc 2 ){AVL 2 + c 2 q 2 ). Using the 
above results for mi(F) and to_i(F) Eq. JJJ reproduces 
exactly the Tkachenko dispersion law 

Let us now discuss the case of a harmonically trapped 
rotating gas, where the non interacting part of the 
Hamiltonian should be replaced by X)fcLi[P 2 /2 TO + 
m(uj 2 _r 2 _+uj 2 z 2 )/2-fl-(rAp)]k E3 with r± = y / x 2 + y 2 . 
A first estimate of the compressional effects on the fre- 
quency of the lowest Tkachenko mode was obtained 
using the dispersion law J^} and setting q — a/ R± , where 
R± is the Thomas-Fermi radial size of the gas |l8j and 
a is a dimensionless parameter characterizing the dis- 
cretization of the normal modes. The value a — 5.45 for 
the lowest azimuthally symmetric mode was extracted 
from the result of Ref. [£| holding in the incompress- 
ible regime where the frequency is linear in q. Since 
in a trapped gas the relevant excitations are affected 
by the compressibility of the gas, it is not obvious that 
the above estimate is enough accurate and it is hence 
important to have more precise calculation of the fre- 
quency of the lowest Tkachenko mode, taking into ac- 
count the finite size, the inhomogeneity as well as the 
compressibility of the gas. In the following we will pro- 
vide such an estimate using the method of sum rules. We 
will make the ansatz F = J2k=i [P( r ±)^z + Q{ r ±)]k for 
the operator exciting the lowest Tkachenko mode. Here 
£ z is the z-direction component of the angular momen- 
tum operator and P and Q are real functions. Following 
the same considerations made in the uniform case, we 
find that the average value of the double commutator 
[Ft, [H, F]] = (h 2 /m) Ek=i( p,i z + Q')l exactly vanishes 
in the hydrodynamic approximation, where £ z — mr A v 
with v = f2 A r at equilibrium, provided the condition 
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is satisfied. As in the uniform case, one can then intro- 
duce the unitary transformation U = e l0F , which pro- 
duces the vortex displacement e(0) = WPe z A r. Ex- 
panding E(9) = (U~ 1 HU) in 9, the mi(F) sum rule is 
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FIG. 1: Lowest Tk achenko freq uency of a trapped Bose gas 
in units of u>o = y/ hQ./ 4mR\ as a function of the angular 
velocity Q. The full line is the sum rule result, while the 
dashed line is the prediction of Eq. (JSJ with q = 5.45 ARx (see 
text). Experimental points are taken from Refs. ITU 1 12 1 II us 
the 3D Thomas- Fermi value for R± . The prediction of Ref. 
corresponds to ujt/loo = 5.45. 



given by the second order term of the expansion which 
can be calculated using the elastic energy change |J3J|. 
One finds the result mi {F) = h 2 (hn/8)(P' 2 rl) . The 
to_i(F 1 ) sum rule is instead extracted from the static hy- 
drodynamic response. The variations of the density and 
of the velocity induced by the static perturbation — XF 
take the form 0] 

Sn = [8fl mfiri +8p+ X(PmrLr 2 ± + Q)]/g , (9) 
Sv = (Sn + XP) e z A r , (10) 

where the changes <5f2 and Sfi are obtained imposing 
the conservation of the number of particles and of an- 
gular momentum. After some straightforward algebra 

one finds: m_i(-F) = m ((P 2 ?!) - (Pr 2 ) 2 / (r 2 )^ /2 + 

Jdr (PmQrl + Q)Sn/2X. 

The above results hold for any choice of the functions P 
and Q satisfying the condition (JSJ . We have written P as 
a polynomial expression of the form P — p s (r±/ R±) s 
and analogously for Q. Using the Thomas- Fermi equilib- 
rium profile ra = {[i— [(w 2 — fl 2 )r^_ + uj z z 2 ]/2} /g, where 
p = m(w 2 _ — Q 2 )R 2 L /2 is the chemical potential, all the 
integrals involved in m\(F) and ni-i(F) are analytical 
and one can write the frequency in the form 
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We have determined the frequency of the lowest mode 
by minimizing the function / with respect to the coef- 
ficients p s at fixed fl/ui±. In practice good convergence 
to the exact value of lot is already ensured for s = 3. 
The results are reported in Fig. ^ When £1 <C w± the 
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FIG. 2: Density change of the lowest Tkachenko mode at 
Q/lo± — 0.9. The inset shows the corresponding displace- 
ment e of the vortex lattice at the angles it/6, n/2, 5ir/6. The 
amplitude of the oscillation corresponds to Sfl/Q = 0.05. 



function / approaches a constant value. When f2 — > lo±_ 
it instead vanishes like 1 — Vt 2 /lo\. Notice, however, 
that if £1 is too close to to± the validity of the present 
calculation, based on the Thomas-Fermi approximation, 
breaks down due to the diluteness of the gas produced 
by the rotation. In Fig. ^ we also show the prediction 
obtained using Eq. @ with q = a/R± and a = 5.45. 
Using the relation mc 2 — /i it is immediate to verify that 
also in this case the dispersion law can be factorized in 
the form Q} with f(x) = (1 - x 2 )/^- 1 ^ - x 2 ) + b^x 2 } 
and a = a 2 = 29.7, b = a 4 /8 = 110.3. The results of the 
sum rule approach instead correspond to a = 31.3 and 
b = 75.8 |20j. The relative difference between the two 
predictions is more and more pronounced as — > lo±_ 
pointing out the inadequacy of the choice q = 5A5/R± 
when R± is large and q becomes small. Fig. ^ also shows 
that the sum rule prediction is systematically closer to 
the experimental values and that the deviations of the 
measured Tkachenko frequencies from the Thomas-Fermi 
values become important only for values of f2 very close 
to where quan tum Hall effects should be taken into 
account HEUEE!- 

In Fig. [21 we report the shape of the density deforma- 
tion © associated with the Tkachenko oscillation. A 
density perturbation of this form, produced by a proper 
change of the trapping potential, should result in a sig- 
nificant excitation of the Tkachenko mode. This density 
change differs from the scaling deformation associated 
with the radial breathing hydrodynamic mode. Actually 
the breathing and the Tkachenko modes are orthogonal 
as confirmed by the fact that the density variation (J5J) 
does not result in any change of the average square ra- 
dius: Jdrr^Sn = 0. In Fig. (inset) we also show 
the amplitude of the corresponding vortex lattice defor- 
mation, obtained from the relation e = Sv/ujt- Similar 



shapes have been obtained in the theoretical calculations 
ofRefs. 0,0 and observed experimentally in Refs. [9ll2l|. 

In conclusion we have shown the occurrence of impor- 
tant compressional features exhibited by the Tkachenko 
oscillations in rotating Bose gases. Sum rules have per- 
mitted to provide accurate estimates of the frequency of 
the lowest mode allowing for a detailed comparison with 
experiments. 

Useful discussions with Gordon Baym, Jean Dalibard 
and Peter Engels are acknowledged. 
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